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We carry out a direct numerical simulation (DNS) study that reveals the effects of polymers on 
statistically steady, forced, homogeneous, isotropic fluid turbulence. We find clear manifestations 
of dissipation-reduction phenomena: On the addition of polymers to the turbulent fluid, we obtain 
a reduction in the energy dissipation rate, a significant modification of the fluid energy spectrum, 
especially in the deep-dissipation range, signatures of the suppression of small-scale structures, 
including a decrease in small-scale vorticity filaments. We also compare our results with recent 
experiments and earlier DNS studies of decaying fluid turbulence with polymer additives. 

PACS numbers: 



I. INTRODUCTION 



The addition of small amounts of polymers to a tur- 
bulent fluid leads to dramatic changes that include mod- 
ifications of the small-scale properties of the flow [J-Q 
and, in wall-bounded flows, the phenomenon of drag re- 
duction [D-[6| , in which polymer additives allow the main- 
tenance of a given flow rate at a lower pressure gradient 
than is required without these additives. Several exper- 
imental, numerical, and analytical studies have investi- 
gated drag reduction (5l4l3j in wall-bounded flows. These 
studies have shown that the addition of polymers mod- 
ifies the turbulence significantly in the region near the 
wall; and this leads to an increase in the mean veloc- 
ity in the bulk. By contrast, there have been only some 
investigations of the effects of polymer additives on ho- 
mogeneous, isotropic turbulence. Examples include re- 
cent experiments jl4l4l6| . which have been designed to 
obtain a high degree of isotropy in the turbulent flow, 
and shell-model studies and direct numerical simulations 
(DNS) 01H013. These studies have shown that 
the addition of polymers to a turbulent flow leads to 
a considerable reduction in small-scale structures; and 
they have also discovered the phenomenon of dissipation 
reduction, namely, a reduction in the energy dissipation 
rate e, in decaying turbulence In this paper we first 

elucidate the phenomenon of dissipation reduction for the 
case of statistically steady, homogeneous, isotropic tur- 
bulence with polymer additives; we then study the small- 
scale properties of such flows. 

We do this by conducting a series of high-resolution 
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DNS studies of the three-dimensional Navier-Stokcs 
equation coupled to an equation for the polymer con- 
formation tensor, which describes the polymer additives 
at the level of the FENE-P @, Q model. Before we give 
the details of our study, it is useful to summarise our 
principal results; these are in two parts. 

The first part contains results from our DNS of forced, 
statistically steady, fluid turbulence, with polymer addi- 
tives, at moderate Reynolds numbers (Re\ ~ 80). The 
forcing is chosen such that the energy injected into the 
fluid remains fixed [liij . both with and without poly- 
mers; this mimics the forcing scheme used in the exper- 
iments of Ref. [ill [la ]. We find that, on the addition 
of polymers, the energy in the statistically steady state 
and the energy-dissipation rate are reduced. This dissi- 
pation reduction increases with an increase in the poly- 
mer concentration c at fixed Weissenberg number We, 
the ratio (see Table [TJ) of the polymer time scale Tp to 
a shearing time scale in the turbulent fluid. The dis- 
sipation reduction also increases with We if we hold c 
fixed. The dissipation reduction seen in our simulations 
should not be confused with the phenomenon of drag 
reduction seen in wall-bounded flows. In the fluid en- 
ergy spectrum we find that the energy content increases 
marginally at small wavevectors on the addition of poly- 
mers, but it decreases for intermediate wavevectors. In 
this part of our study we use 256 3 collocation points and 
attain a moderate Reynolds numbers Re\ ~ 80; but we 
do not resolve the deep-dissipation range. (We consider 
the deep dissipation range in the following paragraph.) 
We also obtain the structural properties of the fluid with 
and without polymers and show that polymers suppress 
large-vorticity and large-strain events; our results here 
are in qualitative agreement with the experiments of 

Ref. mm- 

Furthermore, we find, as in our study of 
decaying turbulence that the polymer extension in- 
creases with an increase in the polymer relaxation time 
rp. We compare our results, e.g., those for the energy 
spectrum, with their counterparts in our earlier study of 
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decaying fluid turbulence with polymer additives Q. In 
such comparisons, we use averages over the statistically 
steady state of our system here; and, for the case of de- 
caying turbulence, we use data obtained at the cascade- 
completion time, at which a plot of the energy dissipation 
rate versus time displays a maximum. 

In the second part of our study we carry out the 
highest-resolution DNS, attempted so far, of forced, sta- 
tistically steady, fluid turbulence with polymer additives; 
we drive the fluid by an external, stochastic force as in 
Ref. H^j. This part of our study has been designed to 
uncover the effects of polymers on the deep-dissipation 
range, so the Reynolds numbers is small (Re\ ~ 16). 
By comparing fluid energy spectra, with and without 
polymers, we find that the polymers suppress the en- 
ergy in the dissipation range but increase it in the deep- 
dissipation range. Finally, we calculate the second-order 
velocity structure function 5*2 (r) directly from the en- 
ergy spectrum via a Fourier transformation; this shows 
that 62(7) with polymers is smaller than 5*2 (r) without 
polymers in this range. 

The remaining part of this paper is organised as fol- 
lows. In Section [TT] we present the equations we use for 
the polymer solution and describe, in subsection III A[ 
the method we use for the numerical integration of these 
equations. Section lTTll is devoted to a discussion of our re- 
sults; as we have mentioned above, these are divided into 
two parts; the first part is contained in subsections IIII Al - 
IHICI and the second in subsection IIII Dl Section |IV] con- 
tains a concluding discussion. 

II. EQUATIONS AND NUMERICAL METHODS 



figurations), I the identity tensor with elements 5 a p, 
f(r P ) = (L 2 - 3)/(L 2 - r 2 ) the FENE-P potential that 
ensures finite extensibility, ?'p = \/Tr{C) and L the 
length and the maximum possible extension, respectively, 
of the polymers, and c = \ij(y + fi) a dimensionless mea- 
sure of the polymer concentration [18j ; c = 0.1 corre- 
sponds, roughly, to lOOppm for polyethylene oxide [f|. 
Table U lists the parameters of our simulations. 
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TABLE I: The cube root TV of the number of collocation 
points, the time step St, the maximum possible polymer ex- 
tension L, the kinematic viscosity v, the polymer-relaxation 
time Tp, and the polymer concentration parameter c for our 
four runs NSP - 256A, NSP - 256B and NSP - 512. We also 
carry out DNS studies of the NS equation with the same nu- 
merical resolutions as in our NSP runs. The Taylor-microscale 
Reynolds number Rex = \/20£ f / \/3vel and the Weissenberg 
number We = tp \fe\Jv are as follows: NSP — 256A and 
NSP - 256B: Rex ~ 80 and NSP - 512: Rex ~ 16; the Kol- 
mogorov dissipation length scale r\ = (v 3 /ef,) 1 ■ For our runs 
NSP - 256A - B, 77 ~ 1.075a;; and for run NSP - 512, 77 ~ 19<5ai, 
where 5x = L/7V is the grid resolution of our simulations. The 
integral length scale l int = (3ir/4) ]T k~ 1 E(k)/(J2 E(k)) and 
T eddy = u rms Ant are as follows: NSP - 256A and NSP - 256B: 
Zmt ~ 1.3 and T od d y =± 1.2; and, for NSP - 512, Z in t ~ 2.05 and 
Teddy ^ 4.0. 



We model a polymeric fluid solution by using the 
three-dimensional, Navier-Stokcs (NS) equations for the 
fluid coupled with the Finitely Extensible Nonlinear 
Elastic-Petcrlin (FENE-P) equation for the polymer ad- 
ditives [|| . The polymer contribution to the fluid is mod- 
elled by an extra stress term in the NS equations. The 
FENE-P equation approximates a polymer molecule by 
a nonlinear dumbbell, which has a single relaxation time 
and an upper bound on the maximum extension. The 
NS and FENE-P (henceforth NSP) equations are 

D t u = ^V 2 u+-^V.[/(r P )C]- Vp + f; (1) 

TP 



D t C = C.(Vu) + (Vu) T .C- /(rp)C 2 . 

Tp 



(2) 



Here u(x, t) is the fluid velocity at point x and time t, 
incomprcssibility is enforced by V.u = 0, D t = dt + u.V, 
v is the kinematic viscosity of the fluid, (i the viscos- 
ity parameter for the solute (FENE-P), Tp the poly- 
mer relaxation time, p the solvent density (set to 1), p 
the pressure, f(x, t) the external force at point x and 
time t, (Vu) T the transpose of (Vu), Cap = (RaRp) 
the elements of the polymer-conformation tensor C (an- 
gular brackets indicate an average over polymer con- 



A. Numerical Methods 

We consider homogeneous, isotropic, turbulence, so we 
use periodic boundary conditions and solve Eq. ([T]) by 
using a pseudospectral method [U [22|. We use N 3 
collocation points in a cubic domain (side L = 27r) . 
We eliminate aliasing errors by the 2/3 rule [U, |22| . 
to obtain reliable data at small length scales; and we 
use a second-order, slaved, Adams-Bashforth scheme for 
time marching. In earlier numerical studies of homo- 
geneous, isotropic turbulence with polymer additives it 
has been shown that sharp gradients are formed dur- 
ing the time evolution of the polymer conformation ten- 
sor; this can lead to dispersion errors [3 [23j]. To 
avoid these dispersion errors, shock-capturing schemes 
have been used to evaluate the polymer-advcction term 
[(u • V)C] in Ref. (23[. In our simulations we have mod- 
ified the Cholesky-decomposition scheme of Ref. [l8| . 
which preserves the symmetric positive definite nature 
of the tensor C. We incorporate the large gradients 
of the polymer conformation tensor by evaluating the 
polymer-advection term [(u • V)£l via the Kurganov- 
Tadmor shock-capturing scheme 24J . For the derivatives 
on the right-hand side of Eq. © we use an explicit, 
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fourth-order, central-finitc-diffcrcncc scheme in space; 
and the temporal evolution is carried out by using an 
Adams-Bashforth scheme. The numerical error in rp 
must be controlled by choosing a small time step St, 
otherwise rp can become larger than L, which leads to 
a numerical instability; this time step is much smaller 
than what is necessary for a pseudospectral DNS of the 
NS equation alone. Table Q] lists the parameters we use. 
We preserve the symmetric-positive-definite (SPD) na- 
ture of C at all times by using [l8| the following Cholesky- 
decomposition scheme: If we define 



J = /(rp)C, 

Eq. ([2]) becomes 

D t J = J.(Vu) + (Vuf.J - s(J - 1) 
where 

L 2 - 3 + f 
S tpL 2 ' 

d/(L 2 -3)-(L 2 -3+j 2 )(f 



■qJ, 



(3) 
(4) 



3) 



(r P L 2 (L 2 -3)) 
Tr {J), and 

Tr[J.(Vu) + {Vuf.J}. 



C and hence J are SPD matrices; we can, therefore, write 
J = CC T , where £ is a lower-triangular matrix with 
elements £ij, such that 1^ = for j > i, and 



J = 



r 2 
' ii 



hxhx l\ X + £ 2 22 hxhx + t2lh2 | ■ (5) 
\lxxhx lixhx+illhl ^1+^2 + ^3 

Equation ((4]) now yields (1 < % < 3 and Tij = diUj) the 
following set of equations: 
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The SPD nature of C is preserved by Eqs. (HJ) if > 0, 
which we enforce explicitly [l8j by considering the evolu- 
tion of \n(£u) instead of la. 

We resolve the sharp gradients in the polymer con- 
formation tensor by discrctizing the polymer advection 
term by using the Kurganov-Tadmor scheme [24| ■ Below 
we show the discretization of the advection term ud x £, 
where u = («, v, w) and I is one of the components of 
the iu/3] the discretization of the other advection terms 
in Eq. ((6]) is similar. 



ud x l 

Hi+l/2,j,k 



i+X/2,j,k 



Hi+l/2,j,k — Hi_i/2J,k 

Sx 

U i+X/2,j,k[^t+X/2,j,k + ^i+X/2,j,k\ 



a i+l/2,j,k[^t+X 



/2,j,k j+l/2,j,feJ 



Sx 

£i+X,j,k T -^■{d x £)i+X/2±X/2,j,k 



a i+X/2,j,k = \ u i+l/2,j,k\ ; 



(7) 



where i, j, k = 0, ... (N — 1) denote the grid points and 
Sx = 5y = Sz is the grid spacing along the three direc- 
tions. 

We use the following initial conditions (superscript 
0): C TOIl (x) = S mn for all x; and u m (k) = 
P„ m (k)v°(k)exp(t6»„(k)), with m,n = x,y,z, P mn = 
(Smn — k m k n /k 2 ) the transverse projection operator, k 
the wave vector with components k m = (-N/2, —N/2 + 
1, . . . ,N/2) and magnitude k = |k|, 6*„(k) random num- 
bers distributed uniformly between and 2tt, and w°(k) 
chosen such that the initial kinetic-energy spectrum is 
E°(k) = k 4 exp(— 2.0fc 2 ). This initial condition corre- 
sponds to a state in which the fluid energy is concen- 
trated, to begin with, at small k (large length scales); 
and the polymers are in a coiled state. Our simula- 
tions are run for 45T e ddy and a statistically steady state 
is reached in roughly 10T c dd y , where the integral-scale, 
eddy-turnover time T e ddy = w rms /Zmt, with u rms the root- 
mean-square velocity and ii n t = ^lkk~ 1 E(k)/^ k E(k) 
the integral length scale. Along with our runs NSP — 256A 
and NSP — 256B we also carry out pure-fluid, NS sim- 
ulations till a statistically steady state is reached; this 
takes about 10 — 15T ddy Once this pure-fluid simula- 
tion reaches a statistically steady state, we add polymers 
to the fluid at 27T e ddy; i-e., beyond this time we solve 
the coupled NSP equations [1] and by using the meth- 
ods given above. We then allow 5 — 6T c dd y to elapse, so 
that transients die down, and then we collect data for 
fluid and polymer statistics for another 25T c ddy for our 
runs NSP - 256A and NSP - 256B. 



III. RESULTS 

We now present the results that we have obtained 
from our DNS. In addition to u(x, t), its Fourier 
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transform u^(t), and C(x, t), we monitor the vor- 
ticity ui = V x u, the kinetic-energy spectrum 
E(k,t) = Efe-i/2<fe'<fe+i/2 l u k'WI ; thc total kinetic 
energy £{t) = YlkE(k,t), the energy-dissipation-rate 
e v (t) = v^ k k 2 E{k,t), the probability distribution of 
scaled polymer extensions P(r'p/L 2 ), the PDF of the 
strain and the modulus of thc vorticity, and the eigen- 
values of the strain tensor. For notational convenience, 
we do not display the dependence on c explicitly. In sub- 
section |nT3] we present the time evolution of E and e v 
and provide evidence for dissipation reduction by poly- 
mer additives. This is followed by subsections IIIIBI and 
IIII CI that deal, respectively, with the effects of polymers 
on fluid energy spectra and small-scale structures in tur- 
bulent flows. In subsection IIII Dl we examine the modi- 
fication, by polymer additives, of fluid-energy spectra in 
the deep dissipation range. 



A. The Energy and its Dissipation Rate 

We first consider the effects of polymer additives on 
the time evolution of the fluid energy E for our runs 
NSP - 256A and NSP - 256B; this is shown in Fig.Q] The 
polymers are added to the fluid at t = 27T ddy The addi- 
tion of polymers leads to a new statistically steady state; 
specifically, we find for We = 3.5 and We = 7.1, that 
the average energy of the fluid with polymers is reduced 
in comparison to the average energy of the fluid with- 
out polymers. By using Eq. (fTJ, we obtain the following 



tives: 



dE 
dt 



Cpoly 



H~ ^poly T ^inj j 
-j/I J u - V 2 U, 

0{^/ u - v[/(rp)c] } 



(8) 



-/f-u 



(9) 



In the statistically steady state ^ = and the energy in- 
jected is balanced by the fluid dissipation rate e„ and thc 
polymer dissipation e po iy Our simulations are designed 
to keep the energy injection fixed. Therefore, we can 
determine how the dissipation gets distributed between 
the fluid and polymer subsystems in forced, statistically 
steady turbulence. 

Before we present our results for the kinetic-energy dis- 
sipation rate, we first calculate the second order-structure 
function S^fV) via the following exact relation (25^ . 



S 2 (r) = r 

Jo 



POO 


^ sm(kr) 


la 


kr 



E{k)dk. 



(10) 



In Fig. [2] we give a log-log plot of £2(7)1 compensated by 
(r/L) -2 , as a function of r/L. We find that, for small 
r, Sz(r) ~ r 2 , which implies that our DNS resolves the 
analytic range, which follows from a Taylor expansion, 
of S2 (r) [2(| ; this guarantees that energy-dissipation rate 
has been calculated accurately. In Fig. [3] we present plots 
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FIG. 1: (Color online) A plot of the fluid energy E versus the 
dimensionless time t/T eddy (runs NSP-256A and NSP-256B) for 
Weissenberg numbers We = 3.5 (blue circles) and We = 7.1 
(black dashed line) . The corresponding plot for the pure-fluid 
case is also shown for comparison (red line). The polymers 
are added to the fluid at t = 27T 8 ddy- 
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FIG. 2: (Color online) Log-log (base 10) plots of the second- 
order structure function 5*2 (r), compensated by (r/L) -2 , ver- 
sus r/L, for our run NSP-256B (blue square) and for the pure- 
fluid case (red circle). The regions in which the horizonal 
black lines overlap with the points indicate the r 2 scaling 
ranges. 



energy-balance equation for thc fluid with polymer addi- of e v {t) versus i/T^ddy for We = 3.5 and We = 7.1 with 



5 



the polymer concentration c = 0.1. We find that the 
average value of e v decreases as we increase We. This 
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FIG. 3: (Color online) Plots of the energy dissipation rate e„ 
versus t/T cddy (runs NSP-256A and NSP-256B) for Weissenberg 
numbers We = 3.5 (blue circles) and We = 7.1(black dashed 
line). The corresponding plot for the pure fluid case is also 
shown for comparison (red line). The polymers are added to 
the fluid at t = 27T cddy . 



suggests the following natural definition of the percentage 
dissipation reduction for forced, homogeneous, isotropic 
turbulence: 

DR = ^ Mzifj^ x 100%; (11) 

here (and henceforth) the superscripts f and p stand, re- 
spectively, for the fluid without and with polymers; and 
the angular brackets denote an average over the statis- 
tically steady state. Percentage dissipation reduction, 
DR, rises with We; this indicates that ef, increases with 
We. [HI Thus, in contrast to the trend we observed 
in our decaying-turbulence DNS H, DR increases with 
We: For We = 3.5, DR ~ 30% and, for We = 7.1, 
DR ~ 50%. Our interpretation is that this increase of 
DR with We arises because the polymer extensions and, 
therefore, the polymer stresses are much stronger in our 
forced-turbulence DNS than in our decaying-turbulence 
DNS (at least for the Reynolds numbers that we achieved 
in Ref. (3). In Fig. H 

we show the cumulative PDF of 
the scaled polymer extension; this shows clearly that the 
extension of the polymers increases with We. In general, 
the calculation of PDFs from numerical data is plagued 
by errors originating from the binning of the data to make 
histograms. Here instead we have used the rank-order 
method to calculate the corresponding cumulative PDF 
which is free of binning errors [28j . 
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rl/L 2 

FIG. 4: (Color online) Log-log (base 10) plots of the cumu- 
lative PDF P c {rl/L 2 ) versus the scaled polymer extension 
r-p/i 2 for We = 3.5 (blue dashed line for run NSP-256A) and 
We = 7.1 (full black line for run NSP-256B). Note that as We 
increases so does the extension of the polymers. These plots 
are obtained from polymer configurations at t = 60T c( jdy 



B. Energy spectra 

In this subsection we study fluid-energy spectra E p (k), 
in the presence of polymer additives, for two different 
values of the Weissenberg number We and fixed polymer 
concentration c = 0.1 (Fig. [SJ. We find that the energy 
content at intermediate wave-vectors decreases with an 
increase in We. At small wave- vector magnitudes k, we 
observe a small increase in the spectrum on the addition 
of the polymers, but this increase is within our numer- 
ical, two-standard-deviation error bars. Because of the 
moderate resolution of our simulations we are not able 
to resolve the dissipation range fully in these simulations. 
We address this issue by conducting high-resolution, low- 
Reynolds-numbcr simulations in Sec. 1111 Dl 



C. Small-scale structures 

We now investigate how polymers affect small-scale 
structures in homogeneous, isotropic, fluid turbulence; 
and we make specific comparisons with experiments (l5l 
ITtjj . We begin by plotting the PDFs of the modulus of 
the vorticity |a;| and the local energy dissipation rate 
eioc = Y^i j(di u j + djUi) 2 /2 in Figs. [51 We find that the 
addition of polymers reduces regions of high vorticity and 
high dissipation [Figs. ((B}]. Furthermore, we find that, 
on normalising |o>| or ei oc by their respective standard de- 
viations, the PDFs of these normalised quantities for the 
fluid with and without polymers collapse onto each other 
(within our numerical error bars) as shown in Figs. [7] 
Our results for these PDFs are in qualitative agreement 
with the results of Refs. QUEH (see Fig. 2 of Ref. Q1] and 
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iog 10 (k) 

FIG. 5: (Color online) Log-log plots (base 10) of the energy 
spectra E p (k) versus k (runs NSP-256A and NSP-256B) for 
c = 0.1 and We = 3.5 (blue squares) or We = 7.1 (black 
stars); we give two-standard-deviation error bars. The corre- 
sponding pure-fluid spectrum E (k) (red circles) is shown for 
comparison. 



Fig. 3 of Ref. [T6j ) . Earlier high-resolution, large- Re\, 
DNS studies of homogeneous, isotropic fluid turbulence 
without polymer additives (see, e.g., Ref. p9l [301] and ref- 
erences therein) have established that iso-|o;| surfaces are 
filamentary for large values of \uj\. In Fig.[5]we show how 
such iso-|u;| surfaces change on the addition of polymers 
(c = 0.1; We = 3.5 or 7.1). In particular, the addition 
of polymers suppresses a significant fraction of these fil- 
aments (compare the top and middle panels of Fig. [8]); 
and this suppression becomes stronger as We increases 
(middle and bottom panels of Fig. [8]). In addition to 
suppressing events which contribute to large fluctuations 
in the vorticity, the addition of polymers also affects the 
statistics of the eigenvalues of the rate-of-strain matrix 
(Sij = (diUj + djUi)/\/2), namely, A„, with n = 1,2,3. 
They provide a measure of the local stretching and com- 
pression of the fluid. In our study, these eigenvalues are 
arranged in decreasing order, i.e., Ai > A2 > A3. In- 
comprcssibility implies that ^\ A^ = 0; therefore, for an 
incompressible fluid, one of the eigenvalues (Ai) must be 
positive and one (A3) negative. The intermediate eigen- 
value A2 can either be positive or negative. In Figs. (|9]) 
and (fTU|) we plot the PDFs of these eigenvalues. The 
tails of these PDFs shrink on the addition of polymers. 
This indicates that the addition of the polymers leads to 
a substantial decrease in the regions where there is large 
strain, a result that is in qualitative agreement with the 
experiments of Ref. [H (see Fig. 3(b) of Ref. [H). Evi- 
dence for the suppression of small-scale structures on the 
addition of polymers can also be obtained by examin- 
ing the attendant change in the topological properties of 
a three-dimensional, turbulent flow. For incompressible, 
ideal fluids in three dimensions there are two topolog- 
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FIG. 6: (Color online) Semilog plots (base 10) of the PDFs 
P(|tt>]) versus |u>| (top panel) and P(ei oc ) versus ei oc (bottom 
panel), for our run NSP-256B, with [c = 0.1, We = 7.1 (blue, 
dashed line)] and without [c = (full, red line)] polymer 
additives. 

ical invariants: Q = -Tr(A 2 )/2 and R = -Tr(A 3 )/3, 
where A is the velocity-gradient tensor Vu |3j|. Topo- 
logical properties of such a flow can be classified [3l], HH 
by a Q — R plot, which is a contour plot of the joint 
probability distribution function (PDF) of Q and R. In 
Fig. HU we give Q — R plots from our DNS studies with 
and without polymers; although the qualitative shape of 
these joint PDFs remains the same, the regions of large 
R and Q are dramatically reduced on the additions of 
polymers; this is yet another indicator of the suppression 
of small-scale structures. 



D. Effects of polymer additives on 
deep-dissipation-range spectra 

In the previous subsections we have studied the effects 
of polymer additives on the structural properties of a 
turbulent fluid at moderate Reynolds numbers. We now 
investigate the effects of polymer additives on the deep- 
dissipation range. To uncover such deep-dissipation- 
range effects, we conduct a very high-resolution, but low- 
Re x (= 16) DNS study (NSP - 512). The parameters 
used in our run NSP — 512 are given in Table fl] The 
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FIG. 7: (Color online) Semilog plots (base 10) of the scaled 
PDFs P(|w|/er) versus |«|/<t u (top panel) and P(e\ oc /a t ) ver- 
sus eioc/ct (bottom panel), where cr„ and a t are the standard 
deviations for \oj\ and ei oc , respectively, for our run NSP-256B, 
with [c = 0.1, We = 7.1 (dashed line)] and without [c = 
(line)] polymer additives. These plots are normalized such 
that the area under each curve is unity. 



fluid is driven by using the stochastic-forcing scheme of 
Ref. [2(| • In Fig. [12] we plot fluid-energy spectra with 
and without polymer additives. The general behavior of 
these energy spectra is similar to that in our decaying- 
turbulence study Q. We find that, on the addition 
of polymers, the energy content at intermediate wave- 
vectors decreases, whereas the energy content at large 
wave- vectors increases significantly. We have checked ex- 
plicitly that this increase in the energy spectrum in the 
deep-dissipation range is not an artifact of aliasing er- 
rors: Note first that this increase starts at wave vectors 
whose magnitude is considerably lower than the dealias- 
ing cutoff fcmax in our DNS; furthermore, the enstrophy 
spectrum k 2 E p (k), which we plot versus k in Fig.[13l de- 
cays at large k; this indicates that the dissipation range 
has been resolved adequately in our DNS. 

For homogeneous, isotropic turbulence, the relation- 
ship between the second-order structure function and the 
energy spectrum is given in Eq. (fT0|) ]25j. Using this re- 
lationship and the data for the energy spectrum shown 
in Fig. I12[ we have obtained the second-order structure 
function S^r) f° r our run NSP — 512. We find that the 





FIG. 8: (Color online) Constant-|u>| isosurfaces for = + 
2<Tuj at t ~ 60T o ddy without (top panel) and with polymers 
[middle panel We = 3.5 (NSP-256A) and bottom panel We = 
7.1 (run NSP-256B)]; |u> is the mean and the standard 
deviation of 



addition of polymers leads to a decrease in the magnitude 
of S2(r). Our plots for S-i(r) are similar to those found 
in the experiments of Ref. [14J . In our simulations we are 
able to reach much smaller values of r/rj than has been 



14 



possible in experimental studies on these systems 
however, we have not resolved the inertial range very well 
in these runs. Note that the spectra E p (k), with poly- 
mers, and E l (k), without polymers, cross each other as 
shown in Fig. 1121 But such a crossing is not observed in 
the corresponding plots of second-order structure func- 
tions (Fig. rTl"]) . This can be understood by noting that 
S2 (r) combines large- and small-fc parts [33[ of the energy 
spectrum. 



IV. CONCLUSIONS 

Wc have presented an extensive numerical study of the 
effects of polymer additives on statistically steady, ho- 
mogeneous, isotropic fluid turbulence. Our study com- 
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FIG. 9: (Color online) Semilog (base 10) plots of the PDF 
P(Ai) versus the first eigenvalue Ai of the strain-rate tensor 
S for the run NSP - 256B, with [We = 7.1 (blue dashed line)] 
and without [c = (full red line)] polymer additives. These 
plots are normalized such that the area under each curve is 
unity. 
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FIG. 10: (Color online) Semilog (base 10) plots of the PDF 
f (A2) versus the second eigenvalue A2 of the strain-rate ten- 
sor S for the run NSP - 256B, with [We = 7.1 (blue dashed 
line)] and without [c = (full red line)] polymer additives. 
These plots are normalised such that the area under each 
curve is unity. 



plements, and extends considerably, our earlier work [3|. 
Furthermore, our results compare favorably with several 
recent experiments. 

Our first set of results show that the average viscous 
energy dissipation rate decreases on the addition of poly- 
mers. This allows us to extend the definition of dissipa- 
tion reduction, introduced in Ref. Q, to the case of sta- 
tistically steady, homogeneous, isotropic, fluid turbulence 
with polymers. We find that this dissipation reduction 
increases with an increase in the Weissenberg number 
We, at fixed polymer concentration c. We obtain PDFs 
of the modulus of the vorticity, of the eigenvalues A„ of 
the rate-of-strain tensor S, and Q — R plots; we find that 
these are in qualitative agreement with the experiments 
of Refs. [IE El. 

Our second set of results deal with a high-resolution 
DNS that we have carried out to elucidate the deep- 
dissipation-range forms of (a) energy spectra and (b) the 
related second-order velocity structure functions. We 
find that this deep-dissipation-range behavior is akin 
to that in our earlier DNS of decaying, homogeneous, 
isotropic, fluid turbulence with polymers [I|. Further- 
more, the results we obtain for the scaled, second-order, 
velocity structure S2 (r) yield trends that are in qualita- 
tive agreement with the experiments of Ref. [H| . 

We hope that the comprehensive study that we have 
presented here will stimulate further detailed experimen- 
tal studies of the statistical properties of homogeneous, 
isotropic fluid turbulence with polymer additives. 
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